Nonlinear surface impurity in a semi-infinite 2D square lattice 
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We examine the formation of localized states on a generalized nonlinear impurity located at, or 
near the surface of a semi-infinite 2D square lattice. Using the formalism of lattice Green functions, 
we obtain in closed form the number of bound states as well as their energies and probability profiles, 
for different nonlinearity parameter values and nonlinearity exponents, at different distances from 
the surface. We specialize to two cases: impurity close to an "edge" and impurity close to a "corner" . 
We find that, unlike the case of a ID semi-infinite lattice, in 2D, the presence of the surface helps 
the formation of a localized state. 

PACS numbers: 71.55.-i, 73.20.Hb, 03.65. Ge, 42.65.Tg 



An interesting recent development for extended, non- 
linear systems with discrete translational invariance is 
the the concept of "breather" or "intrinsic localized 
mode" , whose existence is the result of a careful balance 
between nonlinearity and discretenesai. These excita- 
tions are thought of as generic to a wide range of dif- 
ferent physical systems, including Josephson junctions^, 
biopolymers^, Bose-Einstein condensates in a magneto- 
optical trap^ and arrays of nonlinear optical waveguides^, 
among others. In nonlinear optics, these excitations are 
known as "discrete solitons" (DS) due to their ability to 
move in a more or less robust manner, when endowed 
with momentum (beam angle). In fact, many theoreti- 
cal predictions made for DS have now been experimen- 
tally verified in optics, causing a surge of activity in this 
field. It is believed that an understanding on the cre- 
ation and propagation of DS under different conditions, 
might have a substantial impact on future telecomunica- 
tion/computing systems. 

When looking for discrete solitons, one notes that in 
the limit of high nonlinearity or high power, the effective 
nonlinearity is concentrated in a few "sites" only and, 
therefore, it makes sense to make the approximation of 
replacing the whole nonlinear system for a simpler one, 
consisting of a discrete linear lattice with a a small non- 
linear cluster, or even a single site embedded in it. The 
simplified system is oftentimes amenable to exact mathe- 
matical treatment, leading to closed-form expressions for 
the relevant energies and nonlinearity parameters, as well 
as providing a bound state spatial profile for the relevant 
amplitudes, be these electronic or optical. This high- 
nonlinearity localized state provides a very good starting 
point when looking for discrete solitons in a more general, 
less restrictive context 42^. 

On the other hand, given the practical need to scale 
down the components of any all-optical system, such as 
waveguide arrays, it becomes important to understand 
how the presence of some realistic effects such as bound- 
aries or surfaces affect the creation and propagation char- 
acteristics of these DS. Discrete surface solitons at the 
edge of a one-dimensional (ID) waveguide array has been 
predicted^ and experimentally observed^. It has been 
shown that the presence of nonlinearity can stabilize the 
surface modes in discrete systems, and give rise to differ- 



FIG. 1: (color online) Nonlinear impurity placed near the 
"side" (a) and near the "corner" (b) of a semi-infinite square 
lattice. Along the dashed lines the amplitude is strictly zero. 



ent types of states localized at or near a ID surface, in a 
vibrational^ or optical context^. 

In this work, we consider surface effects for a sim- 
ple two-dimensional (2D) system consisting of a nonlin- 
ear impurity placed near the boundary of a semi-infinite 
square lattice(Fig.l), and examine the conditions for the 
existence of bound state(s), and compare them to the 
results obtained for ID case. 

The stationary modes of a D-dimensional discrete lat- 
tice in the presence of a single nonlinear impurity located 
at d are obtained from the stationary-state discrete non- 
linear Schrodingcr (DNLS) equation 

-«/3C„ + y ^ C„, + 5ndX|Cnr C„ = (1) 

71. n 

where n is a site of a _D-dimcnsional lattice, V is the 
transfer matrix element, x is the nonlinearity parameter 
and a is the nonlinearity exponent. The sum in ^ is 
usually restricted to nearest-neighbors, but other cases 
have also been considered-^. In the conventional DNLS 
case, a = 2 and x is proportional to the square of the 
electron-phonon coupling at site n, while (3 is the eigenen- 
ergy. In nonlinear optics Eq. QJ describes the transversal 
dynamics of an optical field in an array of weakly coupled 
linear waveguides, in the presence of a single, nonlinear 
(Kerr) waveguide. There, a = 2, C„ is the normalized 



amplitude of the field in the nth waveguide, while V is the 
coupling among waveguides and x is the effective nonlin- 
earity of the "impurity" waveguide proportional to the 
nonlinear Kerr coefficient. Also in this case, (3 in Eq.QJ 
must be understood as the propagation constant for the 
allowed optical modes along the longitudinal coordinate 
of the array. Hereafter, for the sake of definiteness, we 
will work in a condensed matter context, but the results 
obtained can be applied to nonlinear optics, where ap- 
propriate. 



I. LOCALIZED STATES NEAR THE SURFACE 
OF A 2D SQUARE LATTICE 

Let us examine the existence of bound states around 
a single generalized nonlinear impurity located near the 
surface of a semi-infinite square lattice (Fig. 1(a), (b)). 
We follow in this section the Green function procedure al- 
ready described in previous work a^?'^-^ , so that the reader 
already familiar with this formalism can skip this sec- 
tion and proceed directly to the next one. We denote 
by d = (dxydy) the position of the impurity. By nor- 
malizing all energies to the half bandwidth of the infi- 
nite chain case (4y), the dimensionless Green function 
G = l/{z — H) can be formally expanded asi^ 



G = G(") + G(")i/iG(") + G(")iJiG(°)iJiG("' 



(2) 



where G^^' is the unperturbed (x = 0) Green function 
of the semi-infinite lattice and Hi = 7|Gd|"|d)(d|, with 
7 = x/4y. Series Q can be resumed to all orders to 
yield 



|ari(") /^(O) 



l-7|Cd|"G|j':i 



(3) 



where iJmn = (m|G|n). The energy of the bound state(s) 
is obtained form the poles of Gmn, i-e-, by solving 



l = 7|CdrG|j"j(zb), 



(4) 



while the bound state amplitudes Gn are obtained form 
the residues of Gmn aX z — z^,: 



|Cd| 



G'-^l{z,)G'^!i{zb) 



(5) 



Inserting this back into the bound state energy equation 
leads to a nonlinear equation for the eigenenergies: 



G 



(0) 
dd 



a+1 



(^6) 



1 



^ [~GT2{zb)r'^' 



(6) 



.(0) 



The unperturbed Green function Gmn for the semi- 
infinite lattice, can be calculated by a judicious appli- 
cation of the method of mirror images, as we will show 
in the next two sections. 



II. IMPURITY CLOSE TO AN "EDGE" 

We start by placing the impurity near the edge of the 
lattice as depicted on Fig. 1(a). In order to simplify mat- 
ters, we take d = (0,d). Since there is no lattice below 
(0, 0), Gmn should vanish identically along the sites lying 
on the dashed line in Fig. 1(b). This implies. 
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where j is a unit vector in the y-direction and where G. 
refers to the Green function of the infinite 2D square 
lattice. Now, using the translation invariance property 



Gl 



G_, 



and the symmetry G 
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mn 



G^, we have 



Gaa = G'^Q and GJ^_d-2j = '^o:2d+2j o^' using a simpli- 
fied notation, 



G 



(0) 
dd 



G(z;0,0)-G(z;0,2d + 2) 



(8) 



where G(z; m, n) refers to the Green function for an infi- 
nite square lattice 



G(z; m, n) 



— r rlfh r dd> cos(TO0i)cos(n02) 

TT^Jo Jo '^';^-(l/2)(cos(</>i) + cos(</>2)) 

(9) 

(see, for instance, ref>i^). We note that IjHJl is identi- 
cally zero at d = —1. The computation of G{z;0,d) 
and G{z; 0, 2d + 2) can be achieved by using some re- 
currence relationsi^ by means of which, an arbitrary 
Green function G{z;m,n) can be expressed in terms 
of two Green functions only: G(z;0,0) and G(z;l,l), 
where G(z;0,0) = i2/7rz)K[l/z^] and G(z;l,l) = 
{2/ttz){ (2z2 - l)K[l/z^] - 2z'^E[l/z^] ), where K[x] is 
the complete elliptical integral of the first kind: K[x] = 

Jq [I — xsm{(j))'^]^^/^d(f>, and E[x] is the complete el- 
liptical integral of the second kind: E[x] = C [1 — 
xs\\\{(t>)^Y''^d(f). In this way, we have obtained a num- 
ber of non-diagonal Green functions in explicit form (see 
Appendix I). In particular, we have obtained G(z;0, 2), 
G(z;0,4) and G(z;0,6) and G(z;0,6) in closed form, 
needed in Eq.©. We finally insert (|SJ| into the RHS 
of the eigenvalue equation, Eq.JHJ, and solve for zi, nu- 
merically. However, the most important features can be 
already deduced from the structure of Eq.®: In Fig. 2 
we show the right-hand side of Eq.JHJl, for the important 
case a = 2 (standard DNLS) and for different d values. 
For comparison, the case d — > oo has also been included. 
Since it is the intersection of these curves with the hor- 
izontal line 1/7 what determines the existence of bound 
states, we see that in general, for finite d a minimum 
value of nonlinearity 7 is needed to create a bound state. 
An increase past the threshold value creates two bound 
states. One of these tends to approach the band while 
the other departs from the band as 7 is increased. As 
argued before in previous worksiS*ii, the former should 
correspond to an unstable localized state, while the latter 
denotes a stable bound state. 
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FIG. 2; (color online) Impurity close to edge: Right-hand side 
of Eq.Q versus z, for a = 2 and for different distances from 
the "edge". 



In Fig. 3 we show a bound state phase diagram in non- 
linearity strength-nonlinearity exponent space, showing 
the number of bound states, for different positions of the 
impurity. As the impurity is brought more and more in- 
side the lattice, the region in parameter space where two 
bound states are possible increases. In the limit d — > oo, 
the curve where a single bound state is found, touches the 
origin, and coincides with the curve for an infinite square 
lattice computed in a previous worki^, as expected. 

An interesting question now is, how does the critical 
nonlinearity to form a localized state 7c depend on d? 
Such critical nonlinearity value is formally given by the 
inverse of the RHS of Eq.©, evaluated at precisely the 
value of z where the RHS of Eq. ((BJ possesses a maximum. 
In Fig. 4 we show 7^ versus d, for a variety of nonlinearity 
exponents. As d is increased past 3, all curves seem to 
converge pretty quickly to their asymptotic values. 

The situation depicted in Figs. 3 and 4 is qualitatively 
similar to what one encounters when placing a nonlin- 
ear impurity near the edge of a semi-infinite ID latticeiS, 
with a difference, though: In the ID case, for a = 2 
the presence of the surface tended to increase 7c, while 
in our case, the proximity of the "edge" tends to de- 
crease 7c: its presence helps localization of the excita- 
tion. We also observe that, for a given impurity position 
d, the nonlinearity needed to create a bound state in- 
creases with a. This was also observed in ID and the 
explanation is quite general, independent of dimension- 
ality: From Eq.(^ we see that, since |C„| < 1, as a is 
increased, I C„ I" will necessarily decrease, meaning that a 
larger value of 7 will be needed to keep the value of the 
effective impurity strength, 7|Cd|". 
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FIG. 3: (color online) Bound state phase diagrams in nonlin- 
earity parameter-nonlinearity exponent, for impurity placed 
at different distances from the "edge": d = O(top left), d=\ 
(top right), d = 2 (bottom left) and d = 3 (bottom right). On 
the solid curve, there is precisely one bound state. 
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FIG. 4: (color online) Scaled critical nonlinearity for onset of 
bound state as a function of the distance from the nonlinear 
impurity to the "edge" of the lattice, for several nonlinearity 
exponents, ranging from q ~ up to a = 4. 



III. IMPURITY CLOSE TO A "CORNER" 

In this case, the impurity is located near the corner of 
the lattice as depicted on Fig. 1(b). In order to simphfy 
matters, we take d = {d^d)\ i.e., we place the impurity 
along the "diagonal" sites. In this case because the im- 
purity is surrounded by "more surface" than in the pre- 
vious case, one would expect even stronger departures 
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FIG. 5: (color online) Impurity close to corner: Right-hand 
side of Eq.JHJ versus z, for a = 2 and for different distances 
along the diagonal. 
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FIG. 6: (color online) Bound state phase diagrams in nonlin- 
earity parameter-nonlinearity exponent, for impurity placed 
at different (diagonal) distances from the "corner": d = O(top 
left), d = 1 (top right), d = 2 (bottom left) and d = 3 (bottom 
right). 



from the ID results already explored in refsiiStii. Since 
there is no lattice to the left or below (0,0), Gmn should 
vanish identically along the sites lying on the dashed line 
in Fig. 1(a). Thus, 



l^^\d- 



2.267, 7^^^((i = 3) = 3.01 (12) 
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(10) 



^d,(-d^-2,dy) T" '~'d,{-d:^-2,-dy- 

We can recast Ea. (|10|) as 

G|j"J(z) = G{z; 0, 0) -2G(z; 0, 2d + 2)-G{z; 2d + 2, 2d+2) 

(11) 
where G{z;m,n) is given by Eq.Q. On Fig. 5 we show 
the right-hand side of Eq.(|SJ, for the important case 
a = 2 (standard DNLS) and for different d values. For 
comparison, the case d ^ oo has also been included. We 
note an important difference with the case of the previ- 
ous section: As z -^ 1+, the RHS of Eq.JH)) approaches 
a finite, non-zero value. This implies the following: An 
increase past a minimum value of nonlinearity 7c cre- 
ates two bound states. One of these tend to depart from 
the band while the other approaches the band as 7 is 
increased. The former state is stable while the latter 
is unstable and, in fact ceases to exist altogether when 7 
reaches a second critical value 7c , marked with a dot on 
Fig. 5. Afterwards, there is only a single bound state. The 
value of this second critical nonlinearity can be obtained 
in closed form by taking the limit z — > 1^ in Eq.(|Bl)(see 
Appendix II). In this way, we have obtained: 



li'\d 



0) = 1.236, 7?Hd= 1) = 1.626 



As d increases, this critical nonlinearity parameter in- 
creases rapidly, and tends to diverge for d — > 00, that is, 
for an infinite square lattice, the unstable bound state 
will still be present at arbitrarily large nonlinearity pa- 
rameter values, a well-known faclii^. 

In Fig. 6 we show bound state phase diagrams in non- 
linearity strength-nonlinearity exponent space, showing 
the number of bound states, for different positions of the 
impurity. As the impurity is brought more and more in- 
side the lattice, the region in parameter space where two 
bound states are possible increases. In the limit d — > 00, 
the region comprising one bound state will get more and 
more "squeezed" into the 7 axis and will formally disap- 
pear for a truly infinite square latticei^. 

In Fig. 7 we show 7c versus d, for a variety of non- 
linearity exponents. We note that, as the impurity is 
brought closer and closer to the "corner" , 7c increases 
or decreases, depending on whether a is above or below, 
approximately, two. This feature was also present in the 
previous case (see Fig. 4). However, in this case, the prox- 
imity effect of the corner is much more pronounced. In 
particular, as the impurity is brought closer to the corner, 
the nonlinearity needed to create a bound state decreases 
even more than when the impurity is brought closer to 
the edge. This implies an even greater departure from 
the ID results. As d is increased past 3, all curves seem 
to start converging towards their asymptotic values. 
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FIG. 7: (color online) Scaled critical nonlinearity for onset of 
bound state as a function of the distance from the nonlinear 
impurity to "corner" of the lattice, for several nonlinearity 
exponents, ranging from a ~ up to a = 4. 



IV. CONCLUSION 



We have examined the formation of bound states 
around a general nonlinear impurity located at, or near 
the "edge" or the "corner" of a semi-infinite 2D square 
lattice. By means of the lattice Green functions for- 
malism, we have obtained in closed form the nonlinear 
equation for the bound state energies, from which we 
have obtained bound state phase diagrams in nonlinear- 
ity strength-nonlinearity parameter, for different impu- 
rity positions with respect to the surface. In general, one 
finds that, a minimum value of nonlinearity is needed to 
create a bound state. Up to two bound states are possi- 
ble, although only one of them is always unstable. These 
features have been observed previously for the ID semi- 
infinite system. However, for the standard DNLS case 
(a — 2), some interesting departures from the ID case 
were also found: (i) The increased number of surface sites 
surrounding the impurity when it is close to the "cor- 
ner" seem to obliterate completely the unstable bound 
state, for relatively high nonlinearity values, (ii) As the 
impurity is brought closer to the surface, the nonlinear- 
ity needed to create a localized state decreases, specially 
in the case when the impurity is near a "corner" . This 
suggests that, in a more general context, when consider- 
ing the creation of discrete solitons near the surface of a 
completely nonlinear (Kerr) 2D square lattice, the sur- 
face (edge, corner) of the square lattice would exert an 
attractive potential, instead of the repulsive one observed 
in semi-infinite ID systemscS. This would make the cre- 
ation of discrete solitons easier to accomplish and observe 
near the boundaries of 2D discrete periodic systems. 
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VI. APPENDIX I 

For an infinite square lattice there are a number of 
recursion relations that allow one to express any de- 
sired Green function, in terms of , ultimately two basic 
ones. We state some recursion relations (see for instance, 
MoritaM). For the sake of space, we drop mention of the 
re-scaled frequency z inside the argument of the Green 
functions and define: 



ylEEG(0,0) = —K[l/z'] 
nz 



(13) 



B = G(l,l) = — [ {2z^-l)K[l/z^]-2z^E[l/z^] ] (14) 

TTZ 



and use the relations: 



G{1,0) = zG(0,0)-l 



(15) 



G{m+ l,m+ 1) 



4to 
2m + 1 
2m -1 



{2z^ - l)G(m,m) 



2m + l^("'"^'"'~^) (^^) 
G(m,+ l,m) — 2zG{m,m) — G{m,m — 1) (17) 
G(m + l,n) = 4zG(m,, n) — G(m, — 1, n) 

-G{m, n+1)- G{m, n - 1)(18) 
G(m + 1,0) == 4zG(m,0)-G(m-l,0) 

-2G(m,l) (19) 



Using these relations, one obtains: 

G(0, 1) = Az-1 

G(0,2) ^-2B-4:Z + A{-1 + A; 

G(l,2) = 1- Az + 2Bz 



(20) 



(21) 



(22) 



G(0,3) = -1 ~ 12 Bz-^ 16 z^ + Az (-3 + 16z^) (23) 



G(2, 2) = (l/3)(-A -AB + SBz^) 



(24) 



(25) 



G(0,4) = 

{l/3)A{-5 + 192z'^) 

-{8/3){B + 22Bz^ + 6{z + Az^)) (26) 

G(l, 3) = (4/3)(A + 6z-6Az'^ + B((7/4) + 4 z^)) (27) 



G(2,3) = (l/3)(-3 + Az + 2Bz(-7 + 8z2)) (28) G{0,6) 



G(3,3) 



-35 8 (-l + 2z^) (-A-4:B + 8Bz^ 



15 



(29) 



G(l,4) = l + A8z^ + 8Bz{3 + 2z^) 

+A{9z-A8z^) (30) 

G(3,4) = 1 + (1/15)^(11 z- 32 2^) 

+ (4/15)Bz(29-84z2 + 64z4) (31) 

G(4,4) = (64/105)A(-(71/64) + 6z2-6z4) 
+ (8/105)S(-(ll/32) + (59/16) z^ 
-9z'' + 6zS) (32) 

G(0,5) = (16/3)48 2" (-1 + A z) 

+ (l/48)(-3-65ylz-140Bz) 

+ (1/3)^2 (-27+ 15 AZ-5OS2) (33) 

G(2,4) = (l/15)A(-23+156z^)-(180/15)z 

+ (B/15)(-19-136z2 + 96z4) (34) 

G(l,5) = 8z(3 + 32 2^) + (1/15)^(28 + 504 22-3840 2-*) 
+ (1/15)5(43 + 2512 2^ + 768 2*) (35) 



(1/15) A(-31 - 2308 2^ + 11520 2" + 15360 z^) 
-(2/15)(30 2 (9 + 256 z^ + 256 2") 
+(1/15)5(23 + 3472 22 + 8768 2*) (36) 



VII. APPENDIX II 



For an impurity close to a "corner" , there is a criti- 
(2) 
cal nonlinearity value 7c , beyond which, the unstable 

bound state ceases to exist. It can be computed by tak- 
ing the limit 2(, — > 1+ in Eq.©, with the Green functions 
obtained in section II. In this way, we have obtained: 



l?Kd=2) 
7f)(d^3) 



27 



(4-^) 



64 (37r - 8)3 
85757r2(657r - 208) 



1.236 



1.626 



= 2.267 



1024(544 -1757r)3 
1320684757r2(6479557r - 2037676) 

64(5668760 - 18052657r)3 
60871560757r2(746695957r - 234592192) 

16384(132029312- 420269857r)3 
= 3.01 
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